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ABSTRACT 

In this paper, least squares and linear minimum mean squared error estimators for channel matrix estimation 
in Rician fading multi-antenna systems is analyzed. The design of mean squared error minimizing training 
sequences is investigated in Multiple-Input Multiple-Output Rician fading channels. The problem of Multiple - 
Input Multiple-Output channel estimation has mostly been treated within the context of minimizing the mean 
squared error of the channel estimate subject to various constraints, such as an upper bound on the available 
training energy. A more general framework for the task of training sequence design in Multiple -Input Multiple- 
Output systems is introduced, which can also treat the minimization of channel estimators. It is shown that the 
proposed framework can be used to minimize the training energy budget subject to a quality constraint on the 
mean squared error of the channel estimator, a weighted channel estimation mean squared error and the mean 
squared error of the equalization error due to the use of an equalizer at the receiver (or an appropriate linear 
precoder at the transmitter). 
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I. INTRODUCTION 

An important factor in the performance of multiple 
antenna systems is the accuracy of the channel state 
information (CSI) [1]. CSI is primarily used at the 
receiver side for purposes of coherent or semi 
coherent detection, but it can be also used at the 
transmitter side and adaptive modulation. Since in 
communication systems the maximization of spectral 
efficiency is an objective of interest, the training 
duration and energy should be minimized. Most 
current systems use training signals that are white, 
both spatially and temporally, which is known to be 
a good choice according to several criteria [2, 3] . 

However, in case that some prior knowledge on the 
channel or noise statistics is available, it is possible to 
tailor the training signal and to obtain a significantly 
improved performance. This context aims at 
providing a unified theoretical framework that can 
be used to treat the MIMO training optimization 



problem in various scenarios. Furthermore, it 
provides a different way of looking at the 
aforementioned problem that could be adjusted to a 
wide variety of estimation-related problems in 
communication systems. 

In this paper, TBCE method is studied in the flat 
Rician fading MIMO channels. We investigate The 
least square(LS) and linear minimum mean square 
error (LMMSE) estimators for channel matrix 
estimation in Rician fading multi-antenna systems is 
analyzed, and especially the design of mean square 
error (MSE) minimizing training sequences. 

The rest of this paper is organized as follows. Section 
II introduces the system model. The training 
sequence optimization in the Rician fading MIMO 
channels is investigated in Sections III. Simulation 
results are presented in Section IV. Finally, Section V 
presents concluding remarks. 
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II. SYSTEM MODEL 

Let us consider a MIMO system with t transmitter 
and r receiver antennas. It is assumed the block 
fading model for flat MIMO channels. Each 
transmitted block has N sub-blocks which contain 
training and data symbols as shown in Fig. 1. The 
frame structure is the same for all transmitting 
antennas. Training and data symbols are located in 
the first and end part of the sub-blocks, respectively. 
In practice, the channel is estimated using training 
symbols in the training phase. To estimate the MIMO 
channel in each sub-block, it is required that n P > t 
training signals are transmitted by each transmitter 
antenna. The t x n P complex received signal matrix 
can be expressed as 



where y = vec Y ,z = vec n , X = X ® I and 



Y = HX + n 



(1) 



where X and n are the complex t -vector of 
transmitted sequences on the t transmit antennas and 
r -vector of additive receiver noise, respectively. 
In MIMO Rician fading channels with K as Rice 
factor, the r x t matrix of channel, H, is defined in the 
following form: 



H = 



k+l '^ray + 



k + l 



H 



los 



(2) 



The matrix H ray explains the Rayleigh component of 
the channel and the matrix Hios describes the channel 
mean value or the Line of Sight (LOS) component of 
the channel. The elements of the matrix H ra y are i.i.d. 
complex Gaussian random variables with the zero 
mean and the unit variance. 



h = vec H = A l——.H rm) + J——.H. . h has the 

k + l y Hk+l l ° s 

following n * rt correlation matrix: 



R,. =E hh" = 



1 



7i 



l + k 



l + k 



K 



K 



l + k 



(4) 



Using 

h = vec(H) = Jy x 



equations 



+k h ™y + 




1 + 



k Aos and ( 4 )< 



the n x rt co-variance matrix of the Rician fading 
MIMO channel can be written as 



C h =R h -E h E h " =^-I« 

l + k 



(5) 



A) ML Channel Estimator 

In classical estimation, the channel is assumed to be 
unknown deterministic. For linear model of the ML 
estimator which maximizes the joint probability 
distribution function (pdf) of is optimal [5]. In 
classical estimation, the channel is assumed to be 
unknown deterministic. For linear model of the ML 
estimator which maximizes the joint probability 
distribution function (pdf) of is optimal. 



K y:h = 



1 



del C II 



rnp 



exp 



H 



- y-Xh C v y-Xh 



(6) 
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Fig. 1. Frame structure for each transmitting antenna 
in a MIMO channel. 

The MIMO channel model of (1) can be expressed in 
the following vector form: 

y = Xh + n (3) 



For noise vector in (3), the co-variance matrix is 
C =R =E vv H =1 



rnp 



(7) 



B) LMMSE Channel Estimator 

For linear model, the MMSE and LMMSE estimators 
are identical. So, let us obtain a linear estimator that 
minimizes the estimation MSE of H. It can be 
expressed in the following general form: 



H 



= M+ Y-MX A 



LMMSE 



O 



Where A„ has to obtained from 

2 



o 



I - F 

J LMMSE 



H-H 



LMMSE 



F 



(8) 



(9) 
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dJ 

The optimal A o can be found from LMMSE /^a ~ ^ 
and is given by 




A o = X H C H X+a 2 n N R I NP 



X H C 



H 



(10) 



The LMMSE estimator of H can be rewritten in the 
following form using (9) in (7): 



H lMMSE =M+Y-MX X h C h X +( j:N r I np X"C h (11) 



-1 



Note that in the Rayleigh fading channel, M = 0, Ch 
=Rh The performance of this estimator is measured 
by the error matrix e = H — H ^ whose pdf is 
Gaussian with zero mean and the covariance matrix. 

f , V 1 



C=R=E s H s = 



C 1 + 

H 



1 



V 



°nN R 



XX H 



(12) 



J 



Therefore, the error of MMSE estimation can be 
computed as 



/ -F 

J LMMSE ~ u 



H-H 



LMMSE 



2 ( 
= E t£ H £ =t C =t 



v 



<N R 



XX H 



(13) 



III. TRAINING SEQUENCE OPTIMIZATION IN 
MIMO 

Let the training matrix P represent the training 
sequence. This matrix fulfils the total power 

constraint^ P P = p which represents the 

maximal number of spatial channel directions that 
the training can excite. The combined received matrix 
of the training transmission is 
Y = HV+N (14) 

where the combined disturbance matrix N is 
uncorrelated with the channel H. Next, we consider 
the problem of designing the training sequence to P 
optimize the performance of the MMSE estimator in 
(11). The performance measure is the MSE and the 
optimization problem can be formulated as 



H 



mint P T ®I S~ l P T ®I 



-i 



(15) 



Observe that the MSE depends on the training matrix 
P and on the covariance matrices of the channel and 
disturbance statistics, while it is unaffected by the 
mean values. Thus, the training matrix can 
potentially be designed to optimize the performance. 

IV. SIMULATION RESULTS 

In this section, our goal is to compare the 
performance of the LS and LMMSE channel 



estimators in the spatially correlated Rayleigh and 
Rician fading channels, numerically. We assume a 2 x 
2 flat fading MIMO system, i.e, Nr = Nt = 2. The SNR 

P 

is defined as r> , as a performance measure, we 



n 



consider the channel MSE, normalized by the 
average channel energy as 



E 



NMSE = 



H.H 



F 



E H 



(16) 



The performance of the MMSE estimators and the 
training sequence design will be illustrated 
numerically. 

In Fig. 1 we give the normalized MSEs averaged over 
5,000 scenarios with different coupling matrices with 
Nr = Nt = 2. The performance of four different 
estimators with MSE minimizing training matrices 
are compared: the MVU/ML channel 



estimator H = YP H PP 



-1 



The one-sided linear 



estimator in [8] that was incorrectly claimed to be the 
linear MMSE estimator. The MVU/ML estimator is 
unaware of the channel statistics (i.e., non-Bayesian), 
and it is clear from Fig. 1 that this leads to poor 
estimation performance. The two-sided linear 
estimator also performs poorly under the given 
premises, but can provide good performance in 
special cases. 



Results based on theoretical formulas 




6 8 10 12 14 
Total Training Power (dB) 



Fig. 2: The average normalized MSEs of channel 
matrix estimation. The performance of four different 
estimators with MSE minimizing training matrices is 
compared. 

In Fig. 2, the performance of the MMSE estimator is 
shown for a uniform training matrix The MMSE 
estimator with three different training matrices is 
compared with the one-sided linear estimator 



P = 



P 



Nn 



(17) 



MSE minimizing training matrix (achieved 
numerically), and the simple explicit training matrix 
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proposed in Heuristic 1. The one-sided linear 
estimator is given as a reference. We used the 
coupling matrix to describe an environment with two 
small scatterers, two big scatterers, and one large 
cluster. It is clear that the gain of employing an MSE 
minimizing training sequence is substantial, and the 
heuristic approach captures most of this gain 
although uniform training is asymptotically optimal 
at high training power. 



Results based on Monte-Carlo simulations 
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Fig. 3: The normalized MSEs of channel matrix 
estimation as a function of the total training power in 
a system .The MMSE estimator with three different 
training matrices is compared with the one-sided 
linear estimator. 

The average optimal training sequence length (i.e., 
average rank of ) is shown in Fig. 3 for both an MSE 
minimizing training matrix and the training matrix 
proposed in Heuristic 1. The average length is given 
as a function of the total training power and for the 
spatial correlation induced by a G {0.33, 0.66, and 1}. 
The optimal training length decreases and the 
convergence towards full rank become slower. The 
heuristic training approach is clearly overestimating 
the training length. 
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Fig. 4: The average optimal training sequence length 
(smallest length that minimizes the MSE) as a 
function of the total training power p. 

V. CONCLUSIONS 

For training-based estimation of Rician fading MIMO 
channel matrices has been introduced, for the 



purpose of joint analysis under different noise and 
interference conditions. The MMSE estimator was 
analyzed in terms of the MSE minimizing training 
sequence and the optimal training structure. Finally, 
the framework was extended to MMSE estimation of 
the channel, using the same type of training 
sequences as for channel matrix estimation. 
Although the MSE of this estimator can be non- 
convex, the limiting solutions at high and low 
training power were derived. 
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